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In this paper, we study non-interacting bosons in a disordered one-dimensional optical lattice in 
| a harmonic potential. We consider the case of deterministic disorder produced by an Aubry-Andre 

potential. Using exact diagonalization, we investigate both the zero temperature and the finite tem- 
perature properties. We investigate the localization properties by using an entanglement measure. 
£Nj . We find that the extreme sensitivity of the localization properties to the number of lattice sites in 

• finite size closed chains disappear in open chains. This feature continues to be present in the pres- 

' ence of a harmonic confining potential. The disorder is found to strongly reduce the Bose-Einstein 

^""^ , condensation temperature and the condensate fraction in open chains. The low temperature thermal 

depletion rate of the condensate fraction increases considerably with increasing disorder strength. 
We also find that the critical disorder strength required for localization increases with increasing 
strength of the harmonic potential. Further, we find that the low temperature condensate fraction 
undergoes a sharp drop to 0.5 in the localization transition region. The temperature dependence of 
the specific heat is found to be only marginally affected by the disorder. 
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I. INTRODUCTION 



The prediction of the localization of Bloch waves moving through a randomly disordered crystal is one of the 
fundamental results in quantum mechanics 1 . Anderson localization of matter waves is a strong function of the 
system dimensionality. In two and lower dimensions, all the single particle states are localized by non-zero disorder 2 ' 3 . 
Anderson localization of matter waves was directly observed in several experiments recently 4-7 . The localization of all 
the single particles states in one and two dimensions occurs for random disorder. However, if the disorder distribution 
is deterministic, localization occurs only beyond a critical disorder strength. A particularly simple model in which 
such a transition from the extended to the localized states occur is the one dimensional Aubry-Andre model 8 . Among 
the experiments mentioned above, an experiment of particular interest for the purpose of this paper is the studies 5 , 
by the LENS group, of the localization of a Bose-Einstein condensate of non-interacting lattice bosons subjected to 
] the Aubry-Andre potential. Our primary motivation comes from these studies in which a finite one-dimensional (Id) 
^ — , ■ Aubry-Andre model was experimentally realized. There have been extensive theoretical studies of the AA model in 
the past 9 " 15 , and in the recent years 16-25 in response to the experimental developments mentioned above (for recent 
CO . reviews see Refs. 26-28). 

In this paper, our main focus is the finite temperature properties of bosons in finite Id optical lattices with AA 
y—{ , disorder in harmonic confining potentials. The incommensurability parameter we consider is q — (v / 5+l)/2. Before 
£N| ' going to the studies of a bose condensate in the AA potential in a harmonic trap, we look at the sensitivity of 
t-H ' localization properties of a single boson in finite optical lattices with open and closed boundary conditions. We note 
that the open boundary conditions are more appropriate for the experiments involving bosons in optical lattices. 
We will show that the extreme sensitivity of the localization properties of a particle in the AA model with closed 
boundary conditions disappear when open boundary conditions are used. This feature is important when we consider 
the many-boson system in finite lattices since now one may use any lattice size. Though this continues to be the overall 
feature after the application of an additional harmonic confining potential, we find a small derealization tendency 
before the complete localization transition occurs. This feature results from a competition between the influences of 
the harmonic potential and the disorder potential. Further, we go on to study bose condensation in the AA potential. 
We study the effects of increasing disorder on the condensation temperature, condensate fraction, and the specific 
heat of the system. The disorder is found to have significant effects on some of these properties. The disorder strongly 
reduces the condensation temperature and the condensate fraction near the localization transition. We also find that 
the fall in the condensate fraction shifts to larger disorder strengths with increasing harmonic potential strength. 
Since the focus of this paper is the effects of deterministic disorder on lattice bosons in harmonic potentials, it is 
fitting for us to place these studies in the general context of existing work on bosons in disordered environments. In 
the past, the effects of random disorder on interacting 3d continuum bosons have been extensively investigated 29 " 37 . 
It has been found that both the condensation temperature and the condensate fraction decreases with increasing 
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disorder strength. In the case of non-interacting lattice bosons, the condensation temperature was found 38 ' 39 to 
decrease for small filling while showing the opposite trend for large filling. There has been extensive studies 42-45 on 
the effects of disorder on interacting lattice bosons within the frame-work of the disordered Bose-Hubbard model. 
Recent experimental studies 46 of interacting disordered lattice bosons in a harmonic trap find decreasing condensate 
fraction with increasing disorder strength. 

The rest of this paper is organized as follows. The study of the single particle localization properties in the open 
and the closed finite length chains is given in the next section (Sec. II). The study on the influence of the disorder 
strength on the condensation temperature, condensate fraction, and the specific heat is presented in Sec. III. The 
conclusions are given in Sec. IV. 



II. THE EFFECTS OF BOUNDARY CONDITIONS ON LOCALIZATION 




As stated, we consider a Id AA model in the presence of a harmonic potential. The Hamiltonian of the system 
is H = —tJ2 <i j > (c\cj + cjcj^ + ^iV^rii, where t is the kinetic energy gain when a boson hop from site i to 

its nearest neighbor site j in the optical lattice with a lattice constant a, cj is a boson creation operator, V(i] = 
Kx\ + \ Cos(2irqxi) is the potential at site i in which K and A are the strengths of the harmonic and disorder 
potentials, q = (\/5+l)/2 is the incommensurability parameter, and Xi is the locator of the site i. All the energies 
in this paper are measured in units of t. After writing this Hamiltonian in a single particle site basis we 
numerically diagonalizc it to obtain its eigen-energies and eigen-functions for various lattice sizes, disorder strengths, 
and the harmonic potential strengths. The localization properties of the eigen-functions are monitored by calculating 
the Shannon entropy which measure the quantum entanglement 22 . The Shannon entropy is given by: S = 
log2Pi, where pi — \af\ in which the a, is the iih site amplitude of the ground state wave- function \ipa) = X^i a * N)- 
The quantum entanglement is maximum for a fully extended state and is zero for a state localized on a site. 

The results of our calculations for open and closed chains, in the absence of the harmonic potential, are shown in 
Fig. 1. In the closed chain results shown, the entanglement (S) falls abruptly at A = 2 signifying a transition from 
an extended state to a localized state for lattices with 144 and 89 sites, but not for other numbers of lattice sites as 



2 



shown in Fig. 1. This extreme sensitivity disappears in the chain with open boundary conditions, as shown in the 
top panel of the Fig. 1. The effect of the harmonic potential on the localization for various values of the harmonic 
potential strength k = Ka 2 is shown in Fig. 2. It is clear that S continues to maintain the overall feature of the 
k = case. However, we find a small increase in S before the complete localization transition occurs. This can be 
understood by examining the nature of the wave- functions shown in Fig. 3. We have plotted the 





FIG. 2. Entanglement (S) as a function of the AA disorder potential strength (A) for an open chain with 100 lattice sites in 
harmonic potentials of different strengths. The curves shown are for : k = (solid), 0.00001 (long dash), 0.00002 (short dash), 
0.00003 (dot), 0.00004 (long dash-dot), and 0.00005 (short dash-dot). The k is given in units of t. 
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FIG. 4. The condensation temperature (To) as function of A (solid line) for 10000 bosons in a Id lattice of 100 sites. The 
harmonic trap strength k = 0.00001. The dashed line shows the variation of the difference between the energies of the first 
excited state (Ei) and the ground state (Eo). 

square of the amplitude of the ground state wave function as a function of i for 100 sites with k — 0.00005 and for 
different values of A > 1. For this case, the af exhibit two peaks around i = 46 and 54. The heights of these two 
peaks increase with increasing A up to A = 2.15. The entanglement, correspondingly, decreases up to A = 2.15 for 
k = 0.00005. As A increases further, the peak heights at i — 46 and 54 positions decrease, while another two peaks 
appear around z = 33 and 67. For A around 2.24 the peak heights at these four positions become almost equal. On 
further increase in A, the peak heights at i = 46 and 54 decrease while those at i = 33 and 67 increase rapidly. We 
find that the S shows a maximum around 2.24 and then decreases rapidly. It is clear that the small increase in S 
results from a competition between the tendency of the disorder potential to localize the wave function on two sets 
of lattice sites away from the center and tendency of the harmonic potential to bring the particle to the center of the 
trap. This section dealt with the single boson properties. In the next section we consider the many-boson system. 

III. BOSE CONDENSATION IN THE AA POTENTIAL 

In this section we study the bose condensation temperature, the temperature dependence of the condensate fraction 
and the specific heat of harmonically confined lattice bosons in AA potentials. To study these properties for a finite 
system, we first numerically diagonalize the Hamiltonian matrix to obtain the single boson energy levels E^. In terms 
of Ei, the total number of bosons is given by the number equation: N = ^2^L N(Ei), wnere Eo an d E m are the 
lowest and the highest single boson energy levels, and N(Ei) = l/[exp[/3 (Ei — /x)] — 1] in which f3 — 1/ksT with 
ks the Boltzmann constant and T the temperature. The chemical potential and boson populations in the various 
energy levels are calculated using the boson number equation. The specific heat is calculated from the temperature 
derivative of total energy (E to t — YllLo N(Ei)Ei) and the condensation temperature 47 ' 48 (To) is determined by solving 
the number equation after setting Nq = and \x = Eq. The effect of increasing disorder strength on the condensation 
temperature is shown in Fig. 4. The disorder clearly leads to a suppression of To. The maximum reduction occurs in 
the transition region between the extended and the localized states. We also find an interesting connection between 
the variation of the T and the disorder dependence of the energy difference (E\ — Eo) between the first excited state 
and the ground state. The Tq is found to closely track (E\ — Eq). The reduction of (E\ — Eo) leads to an increase in 
the one-boson density of energy states at the bottom of the energy specturm. This leads to the reduction in T . The 
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FIG. 5. The temperature dependence of the condensate fraction for different values of A for 10000 bosons in a lattice of 100 
sites in a harmonic potential of strength k = 0.00001. The results shown are for : A = (solid line), 1.0 (long dash), 2.0 (short 
dash), 2.2 (dot), and 2.5 (dash-dot). The inset shows the low temperature variation for A = 2.5. 

temperature dependence of the ground state occupancy with various disorder strengths is shown in Fig. 5. With 
increasing disorder strength, the thermal depletion rate is found to increase considerably. This is especially significant 
as one moves in to the localization transition region. In the localized limit, the condensate fraction decreases sharply 
with a slight increase in temperature as is clear from the results for A = 2.5. The temperature dependence of the 
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FIG. 6. The temperature dependence of the specific heat for different values of A for 10000 bosons in a lattice of 100 sites in 
a harmonic potential of strength k = 0.00001. The results shown are for : A = (solid line), 1.0 (long dash), 2.0 (short dash), 
2.2 (dot), and 2.5 (dash-dot). 
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specific heat for various values of the disorder strength is shown in Fig. 6. Even though in the range of disorder 
strengths shown the system has undergone the localization transition, the specific heat is found to be only marginally 
affected by the disorder. 
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FIG. 7. The variation of condensate fraction with A for 10000 bosons for different harmonic potential strengths. The lines 
are for 100 sites for k values: 0.00001 (solid), 0.00002 (long-dash), 0.00003 (short-dash), 0.00004 (dot), 0.00005 (long dash-dot), 
and 0.0001 (short dash-dot). The circles are for 300 sites for the same k values. The temperature is ksT = 0.1. 

The effect of increasing disorder strength on the thermal depletion of the condensate fraction was shown in Fig. 5. 
We now go on to a detailed investigation of the disorder induced depletion of the condensate. In Fig. 7, we present the 
variation of the condensate fraction with disorder strength for different harmonic potentials. We find that the increase 
in disorder strength leads to two major effects. The first one is the sharp drop in the condensate fraction to 0.5 in the 
localization transition region and second is the shift of this fall to higher values of disorder strengths with increasing 
strength of the harmonic confinement. The latter result correlates well with the localization properties given in Fig. 
2. This shift is due to larger disorder strengths required for the localization transition when the harmonic potential 
increases. In order to understand the first effect, we examined the energy level differences between the ground state 
and the lower excited energy levels. In Fig. 8, we have shown the energy level difference, as a function of disorder 
strength, between the first excited state and the ground state, and between the second excited state and the ground 
state. When the energy difference between the first excited state and the ground state becomes very very small, both 
the states are almost equally populated, and No/N reduces to about 0.5. Note that the energy difference between the 
first excited and the ground state remains finite, though small, even for large values of A, which is not visible in the 
scale of this figure. Note also that that the E 2 — E shows a minimum in the vicinity of the localization transition 
and then increases fast as the system moves deeper in to the localized regime. 
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FIG. 8. The variation of the energy difference AE between low lying levels with A. The results are for k = 0.00001 : Ei — Eq 
(open circles) and E 2 — E (solid line) and for k = 0.00005 : Ei — E (filled diamonds) and E 2 — E (dash line). 



IV. CONCLUSIONS 

In this paper, we investigated some properties of a single boson as well as many boson systems in finite one- 
dimensional optical lattices with Aubry-Andre disorder in harmonic confinements. We find that the open lattices 
does not show the extreme sensitivity of the localization properties seen in finite closed lattices. This continues to 
be the case in the presence of an overall harmonic confining potential. In the presence of the harmonic potential, 
we find a small derealization tendency before a completely localized state is reached. On the basis of the disorder 
induced changes occurring in the single particle wave-function, we argued that this tendency is due to a competition 
between the effects of the harmonic potential and the disorder potential. For the many-boson system, we studied 
the effects of disorder on the condensation temperature, the condensate fraction, and the specific heat. We find that 
the condensation temperature decreases with increasing disorder strength, the rate of the reduction is maximum in 
the localization transition region. The disorder is found to have a significant effect on the condensate fraction. With 
increasing disorder, the thermal depletion rate is found to increase considerably, especially in the localization transition 
region. In the localized limit, the condensate fraction decreases significantly with a small increase in the temperature. 
We also studied the variation of the low temperature condensate fraction with increasing disorder strength. We find 
that the condensate fraction shows a sharp drop to 0.5 in the localization transition region. We also find that this fall 
shifts to higher disorder strengths with increasing strength of the harmonic confining potential. The effect of disorder 
on the specific heat is found to be a marginal reduction in a certain range of temperature. 
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